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We derive a set of equations for the wavefunction describing the marginal bound 
state of a single DO-brane with a single D4-brane. These are equations determining the 
vacuum of an = 8 abelian gauge theory with a charged hyper multiplet. We then solve 
these equations for the most general possible zero-energy solution using a Taylor series. 
We find that there are an infinite number of such solutions of which only one must be 
normalizable. We explore the structure of a normalizable solution under the assumption 
of an asymptotic expansion. Even the leading terms in the asymptotic series, which should 
reflect the supergravity solution, are unusual. Through the Spin{5) flavor symmetry, the 
modes which are massive at long distance actually influence the leading behavior. Lastly, 
we show that the vacuum equations can quite remarkably be reduced to a single equation 
involving one unknown function. The resulting equation has a surprisingly simple and 
suggestive form. 
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1. Introduction 



A single DO-brane and a single D4-brane form a marginal bound state [Q. The low- 



energy dynamics of a DO-brane in the presence of a D4-brane is described by a quantum 
mechanical Yang-Mills theory with eight supercharges. The theory has a U{1) vector 
multiplet coupled to a charged hypermultiplet 0. With a single D4-brane, there is only 
a Coulomb branch. The same quantum mechanics appears in the problem of counting 
H-monopole ground states in the toroidally compactified heterotic string 0. While the 
structure of vacuum wavefunctions in marginally bound systems is typically very difficult 
to analyze, this particular theory has a number of simplifying features. The aim of this 
paper is to study the vacuum wavefunction of this + 1-dimensional gauge theory with 
eight supercharges. 

Our goal is to gain insight into a number of issues. For example, how do we go about 
uncovering the structure of threshold wavefunctions? It is not even clear how to formulate 
reasonable questions about a system as complex as the quantum mechanics describing 
many DO-branes. Another major issue is how the full quantum mechanics resolves the 
singularity of the moduli space metric. The vector multiplet contains five scalars x^. For 
large r = |x|, the effective action describing the Coulomb branch dynamics should be a 
reasonable description of the physics Q. The metric on the Coulomb branch is protected 
by supersymmetry and takes the form. 



where is the Yang-Mills coupling constant. We can express in terms of the type IIA 
string scale Ms and coupling constant gs, 



This is the only scale in the theory. In the following sections, we set g^ = 1 for simplicity. 
The tube-like metric (1.1) has a singularity at r = which is resolved by the full quantum 
mechanics. Metrics with a similar structure appear in D1-D5 systems. 

In the following section, we present the supercharges and describe the symmetries of 
the problem. The fiavor symmetry is Spin{5) x SU{2)ji, and the unique vacuum is invariant 
under this symmetry |0] . In section three, we derive the general form of a gauge invariant 
and fiavor invariant wavefunction. A general wavefunction depends on 11 functions of two 




(1.1) 



9^ = 9sMl 
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variables, r and y. As above, r is a radial coordinate for the 5 scalars of the vector 
multiplet. The hypermultiplet has 4 scalars qi, which parametrize the massive directions. 
We take y = \q\. We then derive a set of differential equations that any zero energy 
wavefunction must obey. 

In section four, we analyze the implications of these differential equations. We can 
immediately reduce the number of unknown functions from 11 to 7. These 7 functions 
satisfy 14 first order coupled partial differential equations in 2 variables. We point out 
some intriguing features of these equations: in particular, there is an interesting formal 
method of reducing the number of functions and equations in which the harmonic oscillator 
plays a key role. This method for collapsing the system of equations has similarities to 
the technique used to find non-renormalization theorems 0. We, however, proceed in 
a different direction. We solve the 14 equations exactly by Taylor expanding in the y 
variable. This is an expansion in the massive directions. The first interesting point is 
that the differential equations alone do not determine the solution. There are actually an 
infinite number of zero energy solutions. The 'gauge' degrees of freedom are not a finite 
set of parameters as we might have expected, but an entire function of degrees of freedom. 

Essentially, from the perspective of the Taylor series, the problem boils down to solving 
2 ordinary differential equations in 3 unknown functions. All the other terms in the 
wavefunction are determined by these 3 functions of r. The condition that must uniquely 
specify the solution is normalizability. Rather remarkably, this global condition is strong 
enough to fix an entire arbitrary function. This seems to hint that the kind of principle 
that should underly M theory involves global rather than local constraints. In a vague 
sense, this is reminiscent of holography. 

In section five, we turn to the practical problem of determining the normalizable 
solution. It turns out to be difficult to implement the global constraint of normalizability 
in any nice way. Instead, we expand the solution in an asymptotic series. This is akin to 
solving the M theory equations of motion for the geometry of an M5-brane in a derivative 
expansion. We find that even the leading terms in the solution are unusual. These terms 
should match a supergravity analysis. However, the structure of these terms is strongly 
dictated by invariance under the Spin{5) flavor symmetry. Invariance under Spin{5) is 
not a statement about long distance, moduli space physics. It is a statement that requires 
knowledge of both long and short distance physics because the Spin{5) generators act on 
both the massless and massive degrees of freedom. Somewhat contrary to intuition, we 
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find that the vacuum state for the massive degrees of freedom at large r is a sum of 3 
representations of Spin{5): a spherically symmetric 1, a 5 and a 14. 

The leading terms in the bound state wavefunction \1/ have the form, 

* ~ ;:3^"~l^i > +74e~~l^2 +—^e-—\b3 (1.2) 

where the \bi > are constructed from fermions. It seems unlikely that this asymptotic 
form could have been determined from low-energy considerations alone. In this sense, the 
massive degrees of freedom, through the Spin{5) flavor symmetry, are important even at 
arbitrarily long distances. The structure of \1/ in (1.2) really begs for an interpretation 
both in terms of the supergravity solution of the D0-D4 brane 0, and in terms of the 
DLCQ description of an M5-brane via Matrix theory 

We proceed to compute the general form of \E' in an asymptotic expansion. This takes 
us well beyond supergravity. The corrections to the leading terms (1.2) take the form 
of a perturbation series in the coupling constant g^. Each correction depends on an a 
priori unknown constant, and we give a prescription for determining this constant. This 
amounts to summing up the higher derivative corrections to the supergravity solution for 
an M5-brane. Is the asymptotic solution actually convergent as r—>-07 We know of no non- 
perturbative effects in the abelian gauge theory that could be relevant at short distances. 
However, this does not prove that the solution is convergent. It would be interesting to 
sum up a sufficient number of terms in the asymptotic series to see whether the solution 
is well-behaved as r becomes small. This would clarify how the singularity in the metric 
is resolved from the perspective of a derivative expansion. 

We note that finding the bound state wavefunction in an asymptotic series is much 
like trying to understand M theory in a derivative expansion. In section six, we present a 
quite different reduction of our initial 14 vacuum equations, one that perhaps an M theorist 
might use. The result is quite incredible. The entire problem reduces to solving a scalar 
equation of the form, 

(a + B -V + W^u = 0, (1.3) 

where A = d'^ + dy and m is a particular combination of the functions that appear in the 
bound state wavefunction. The vector field B and potential W are rational functions of r 
and y. Equation (1.3) is both simple and highly suggestive. It would be very interesting to 
find the solution to (1.3) either analytically or numerically. Supergravity and the structure 
of the derivative expansion should emerge from the asymptotics of the resulting solution 
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for the bound state wavefunction. It would also be extremely interesting to generalize this 
reduction to non-abelian gauge theories. This might possibly help us understand how to 
define M theory through Matrix theory 0. It does not seem too unlikely to us that the 
ability to 'deprolong' the initial vacuum equations to get (1.3), as described in section six, 
is tied to supersymmetry and invariance theorems 

There are many additional directions to explore. Some of the simplifying features of 
the D0-D4 system remain when we add more hypermultiplets. However, there will now 
be a Higgs branch and the zero energy wavefunctions will spread in an interesting way 
onto the Higgs branch. Turning on S-fields makes the gauge theory on the D4-brane non- 
commutative [0. Certain choices of S-field should change the asymptotic behavior from 
polynomial decay to exponential decay. It would also be interesting to actually match the 
asymptotic structure of the bound state wavefunction with higher derivative corrections 



to the supergravity solution, like those generated by the R terms [11,12,13] 



2. The D0-D4 Quantum Mechanics 

2.1. The vector multiplet supercharge 

The D0-D4 system is obtained by dimensionally reducing = 1 abelian Yang-Mills 
with a single charged hypermultiplet from six dimensions. The symmetry group consists 
of the R-symmetries0 Spin(5) x SU{2)ji ~ Sp{2) x Sp{1)r. The Hamiltonian is invari- 
ant under the symmetry group while the eight real supercharges transform in the (4, 2) 
representation. 

Let us begin with the vector multiplet which contains the five scalars transforming 
in the (5, 1) of the symmetry group. Let be the associated canonical momenta obeying, 

[x^',p'']=iS^''. (2.1) 

Associated to these bosons are eight real fermions Xa where a = 1, . . . , 8 transforming in the 
(4, 2) representation of the symmetry group. These fermions obey the usual quantization 
relation, 

{XaAb}=Sab. (2.2) 

^ The symmetry group, including both gauge and flavor symmetries, is not globally a prod- 
uct. There are discrete identifications. However, for this analysis we only need the Lie algebra 
generators so we can ignore global identifications. 
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Let be hermitian real gamma matrices which obey, 

{^M^^^l = 2(5'^^ (2.3) 

An explicit basis for these gamma matrices along with a discussion of the symmetry group 
action is given in Appendix A. 

To write the vector multiplet supercharge, we introduce an auxiliary field D which 
transforms as (1,3) under the symmetry group. The D-term is independent of x'^. The 
vector multiplet supercharge is given by: 

Ql = W>^)a + DabXb. (2.4) 

The real anti-symmetric matrix D commutes with because the Sp{l)ii and -5^(2) actions 
commute. The D-term must also satisfy, 

DacD^c = -Sab (2.5) 

It is then not hard to check that: 

{Ql,Ql}^6ab{p^ + \D\'}. (2.6) 
Under a symmetry transformation {U,g) e Sp{2) x Sp{l)ii, we note that 

7V^t/7'^p^C/-\ X^UgX, D^gDg-\ (2.7) 

so that 

Qv^UgQ\ (2.8) 

2.2. The hypermultiplet supercharge 

A hypermultiplet contains four real scalars which we can package into a quaternion 
q with components q'^ where i = 1,2,3,4. This field transforms as (1,2) under the sym- 
metry group. We again introduce canonical momenta Pi satisfying the usual commutation 
relations. 

The hypermultiplet is charged under the U{1) gauge symmetry so we need to deter- 
mine how U{1) acts on q. The q'^ parametrize so the symmetry group acting on the 
hypermultiplet must sit inside. 



50(4) ~ Sp{1)l X Sp{1)r. 
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Gauge transformations and Sp{1)r transformations commute. Therefore, the U{1) gauge 
symmetry sits inside Sp{1)l- We choose to generate the gauge symmetry by left multiph- 
cation on q by /. The hermitian generator of gauge transformations on the bosons is then 
given by, 

Gb = Wi2 + W34, (2.9) 

where 

Wij = QiPj - QjPi. (2.10) 
The superpartner to is a real fermion i/ja with a = 1, . . . , 8 satisfying, 

{i^a.i^b} ^ Sab, (2.11) 

and transforming in the (4, 1) representation. In terms of the operators given in Ap- 
pendix A, the free hypermultiplet charge takes the form 

Qa' ^si,i;bPj. (2.12) 

Note that since the implement right multiplication by a quaternion, they commute with 
7^*. This free charge obeys the algebra, 

|<5a^,Q6^} = SabPiPi- 

Invariance of (2.12) under the U{1) gauge symmetry requires that 

Gf = ~sl,i;ai^b (2.13) 

generate gauge transformations on ip. The total generator of the U{1) gauge symmetry is 
then given by, 

G = Gb + Gf = Wi2 + Ws4- '-^l^s^. (2.14) 
The full hypermultiplet supercharge also includes couplings to the vector multiplet, 

Qa = sibi^bPj + il^'s^sXbil^hX^'qj. (2.15) 

The form of the interaction term in (2.15) is fixed up to an overall constant by symmetry. 
The appearing in the interaction term is needed to ensure that is gauge- invariant. 
The charge obeys the algebra: 

[Qt Qt] = Sab I^P^P^ + ' \ O^W^V} + 27aV^G'. (2.16) 

As we expect, the supersymmetry algebra only closes on the Hamiltonian up to gauge 
transformations. 
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2.3. The coupled system 

The full supercharge Q is the given by, 

Q = QVj^Qh^ (2.17) 

where we define the D-term in the following way: using the components and matrices 
s*, we can define a quaternion and its conjugate which act by right multiplication on A: 

qR = gigi + s^q^ + s^q^ + s^q'^ 
qR = - - s\^ - s\\ 

We can then write the D-term in the form, 

Dab = \{q''s^q'')^,. (2.18) 

This D-term obeys (2.5), 

= ^ab 2 {qqf ■ 

The full charge obeys the algebra: 

{Qa, Qb} = Sab + \D\^ + P^P^ + \x\^\q\^ + ...}+ 2-f^,X^G 

-5a,2if + 27>^G'. 



^2.19) 



The omitted terms are bilinears in the fermions whose exact form we will not need. The 
bosonic potential V appearing in (2.19) is given by, 

V = + 1 |g|4. (2.20) 

Since we have coupled a single hypermultiplet to the U{1) vector multiplet, the only fiat 
direction is = and there is no Higgs branch. 



3. Deriving Equations for the Vacuum Wavefunction 

To be consistent with predictions from string duality, there should be a unique vacuum 
wavefunction for this quantum mechanical gauge theory [Q. An index argument proves 
that there is at least one normalizable vacuum wavefunction [|l|. Coupled with a recent 
invariance theorem 0, the index result implies that the ground state is unique. 
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On quantization, the fermions A and t/j act as gamma matrices on a 2 5 6- dimensional 
spinor wavefunction. A priori, the vacuum wavefunction then consists of 256 complex 
functions of the 9 bosonic variables and q"^. However, we can significantly simplify the 
problem by using symmetries. First note that any state in the Hilbert space \s > must be 
gauge-invariant , 

G'|s>=0. (3.1) 

Further, all states can be grouped into representations of the global Sp(2) x Sp{l)ii sym- 
metry group. The Sp{2) is generated by the operators, 

= - ^ i^aXb + i^ai^b) , (3.2) 

where 

X^'^ = xV-a;V- (3.3) 

The three generators of Sp{l)ii correspond to right multiplication by /, J, K and in accord 
with prior notation, we will denote them by s*: 

52 = Wi2 -Ws^ + ^ Xs^X 

5^ = Wi3 + + I Xs^X (3.4) 
§4 = -W23+'- Xs^X. 

The unique ground state ^ must be invariant under the actions of Qa, G^T^'^ and 5*, 

Qa^ = G^ = T'"'<i> = 5'<if = 0. 

These constraints are quite powerful; for example, they allow us to replace a differential 
operator X'^'^ by an algebraic one | 7^^ (A^Ab + i/^ai^b) ■ There are multiple ways to derive 
equations for the vacuum wavefunction. We will describe two approaches which we used 
to derive these equations.@ 

2 Using two different approaches helped enormously in the search for errors. 
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3.1. Radial coordinates 

In the first approach, we can rewrite the supercharge (2.17) in terms of radial coordi- 
nates: 

2 I |2 2 I |2 

r =\x\ , y =\q\ . 

The charge takes the form, 

+ {Ys's^)aX^q^ + {D\)a (3.5) 
r y 

We have lumped the angular derivatives and non-derivative terms into the operator M^. 
What is particularly nice about Ma is that we can replace all the derivative operators by 
bilinears in fermions. As we noted before, Sp{2) invariance allows us to replace X^^^ by 

\ lib i^ah + Ipai^b) ■ 

However, we can also replace q'^Wkj by a bilinear in fermions using Sp{l)ji invariance. 
Using (3.4), we note that: 

qf'Wki = - {g^As^A + q^Xs^X + g^As^A} 

q'Wk2 = ^ {q'Xs'X - q'Xs'X - q'Xs^X} 

^ (3.6) 
q'Wks = {q'Xs'X - q'Xs^X - q'Xs'X} 

qf'WkA = ^ {q^Xs^X - q^Xs^X - q^Xs^X} . 

Therefore, Ma is a completely algebraic operator. 

Using the Sp{2) symmetry, we can then rotate x to the special point where 
and x^ = for n > 1. Likewise, we can rotate q using Sp{1)r to the point q^ ^ and 
= for i > 1. At this point, r = |a;^| and y = \q^\. Since all angular derivatives in 
Qa are replaced by algebraic operators, there is no difiiculty in restricting ^ to this point. 
The question of determining at this point then reduces to finding coupled differential 
equations in two variables. The form of ^ at an arbitrary choice of x and q can then be 
obtained by applying the rotation generators (3.2) and (3.4). 
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3.2. Symmetries and the fermion Hilbert space 

In the second approach which we wiU use for the rest of the paper, we will first solve 
the invariance conditions explicitly for the most general possible invariant wavefunction. 
As we will show, the most general wavefunction depends on 11 functions of r and y. 
We will then derive coupled equations for these functions from the requirement that the 
wavefunction have zero energy. 

The first step is to construct the fermion Hilbert space. We need to complexify our 
real fermions and build a Fock space: 





= Ai 


+ 1X2 


V2dvi 


= V'l 


+ #2, 


\/2du2 


= A3 


— iX^ 


V2dv2 


= ^3 


- #4, 


V^du^ 


= A5 


+ iX6 


V2dv3 


= tpb 


+ ^■06, 


V2du4 


= Ar 


- iXg 


V2dv4 


= V'T 


- lips- 



(3.7) 



It is natural to think of dua and dva as one-forms obeying the relation, 

{dUa, dul} = 5 ah {dVa, dvl) = 5 ah, 

where a = 1, . . . , 4. Wavefunctions in the Hilbert space are then (p, q) forms where p and 
q are the dua and dva degrees, respectively. We choose the Fock vacuum or (0, 0) form to 
satisfy, 

du*a\Q >= dva\0 >= 0. 

Note that the complex conjugate of a {p, g)-form is a (4 — p,4: — g)-form. If the ground 
state is unique then it is bosonic so the form degree must be even. 

With our choice of complexification (3.7), the Sp{2) generators T^^^ acting on forms 
preserve degree. Actually, the generators preserve p and q separately so a (p, q) form is 
mapped to a [p, q) form. The Sp{2) generators naturally split into commuting generators 
for an Sp{2)h acting on bosons and an Sp{2)f acting on fermions. In turn, the Sp{2)f 
splits into an Sp{2) acting on du with generators, 

-| J2'^ab>^aXb, 
a,h 

and an Sp{2) acting on dv with generators, 

a,h 
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We can now employ some group theory to see how the various 128 bosonic forms transform 
under Sp{2)f. Let us start with the (p, 0) forms which appear in the foUowing representa- 
tions: 



p 


Sp{2)f rep. 





1 


1 


4 


2 


50 1 


3 


4 


4 


1 



We wedge the (odd) even {p, 0) forms with the (odd) even (0, q) forms to get the 128 
bosonic forms. The foUowing representations appear from wedging even forms with even 
forms, 

(1)^°© (5)*^ 10 14, 
while from wedging odd forms with odd forms, we find: 

(1)"0(5)^©(1O)^ 

We can immediately discard forms transforming in the 10 representation. A tensor say a^'^ 
transforming in the 10 is antisymmetric in u so contraction with x^x'^ to get a singlet 
of the full Sp{2) gives zero. 

Let us now constrain our Hilbert space further by imposing invariance under Sp{l)ji 
and the gauge symmetry. We can rewrite the generators (3.4) in terms of our complex 
fermions: 

a 

= VFi3 + VF24 — i {du\du2 + du[du2 + dusdu^ + duldul) (^"^^ 
s"^ = Wi4 — W23 + {duidu2 — du\du2 + du^duA — du^^du^) ■ 
Likewise for the gauge symmetry, 

G = Wi2 + W34 - J2 + 2- (3-9) 

a 

Note that the operator Wij has eigenvalues n and —n with corresponding eigenfunctions, 
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What does invariance under (3.8) and (3.9) imply? By taking the sum and difference of G 
and s^, we see that we should restrict to (p, q) forms 

(.12)^(.34)^-^|P,^>, (3.10) 

which can be multiplied by a function of both |-2i2p and |-234p. The remaining two gen- 
erators 5"^ and change the value of p. It is natural to study the complex combinations 
— is'^ and + is^, which raise and lower the value of p: 

s"*" = — is^ = Z12P34: — ZS4P12 — 2i {duidu2 + duzduA 

(3.11) 

s = + is"^ — Z12P34 — Z34P12 — 2i {duldu2 + du^dul) , 

where Pij = Pi — ipj and Pij = Pi + ipj. 

Again the generators of Sp{1)r split into an Sp{l)b acting on bosons and an Sp{l)f 
acting on fermions. It is easy to see how the {p, q')-forms fall into representations of -5^(1) /. 
The three singlets under Sp{2)f^ denoted \0,q >, \2,q >ip, \A,q > transform in the 3 of 
Sp{l)f. For the choice = 0, 4, we can construct one singlet under the full Sp{1)r which 
can be multiplied by an arbitrary function of y. For the case q = 2, we can construct two 
singlets under Sp{l)u by tensoring with either the 1 or the 5 of Sp{2)f^. Let us denote 
the 5 of Sp{2)f^ by |2 >5^- The explicit Sp{l)ji singlets are then given by the forms, 

{^34 + Zi2Z34{duidU2 + dU3dU4) + Zi2dUidU2dU3dU4^ |0 >, 

{Z12Z34 + ^(|2;i2p - \z^4\'^){duidu2 + du^du4) - Zi2z-iiduidu2du^du4] x 
{dvidv2 + dv3dv4)\0 >, 

{Z12Z34 + ^(|^12p - \zs4\'^){dUidU2 + dUsdu4) - Zi2Z34dUidU2dUsdU4]\2 >^_^, 
{^12 ~ Zi2Z34{dUidU2 + dU2dU4) + Z24dUidU2dU3dU4^ dvidv2dv2dv4\0 > . 

The 5 oi Sp{2) denoted |2 >5^ decomposes into five singlets under Sp{l)f. The form 
1 2 can therefore only appear with a function of y. To satisfy the constraint (3.10), 
we must then tensor |2 >5^ with a g = 2 form constructed from dv. The two choices are 
either the 1 or the 5 of Sp{2)f^. This gives three additional possibilities denoted, 

|2,2>i, |2,2>5, |2,2>j4, 

where the subscript denotes the representation under the full Sp{2)f. The construction of 
these forms is described in Appendix B. 
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Lastly, we need to consider the case of odd p. The |1 and |3 >4^ forms combine 
to form a doublet under Sp{l)f. By tensoring with either the |1 >4^ or the |3 >4^ forms, 
we can construct the following four Sp{l)ji invariants: 

{Zl2 - Zs4{dUidU2 + dU3dU4)} |1, 3 >i 
{zi2 — Zs4idUidU2 + dusdu4)} |1,3 >5 
{Z34 + Zi2{dUidU2 + dU3dU4)} |1, 1 >1 
{Z34 + Z\2{duidU2 + dU3dU4)} |1, 1 >5 . 

Again the subscript denotes the representation under the full Sp{2)f. After imposing all 
the invariance constraints, we are therefore left with 11 complex functions fi = fi{r,y) 
appearing in the following way: 

/i|0, >= /i 1^34 + zi2Zs4{duidu2 + dusdu4) + zl2duidu2dusdu4^ |0 > (3-12) 



/2IO, 4 >—f2 {-2i2 — zi2Zs4{duidu2 + du2,du4) + Z34duidu2du2,du4^ 
X dvidv2dv3dv4\Q > 



(3.13) 



/3|2,2>i, /4X^|2,2>g, /5X^x'^|2,2>^^, (3.14) 

(3.15) 



/ell, 3 >= fe{zi2 - Z34{duidu2 + du3du4)} |1,3 >i 

/7II, 1 >= /7{^34 + Zi2{dUidU2 + dusdu4)} |1, 1 >1 



(3.16) 



fsx^\l, 3 >^= fgx^ {zi2 - Z34{duidu2 + du3du4)} |1, 3 >g 

hx^\l, 1 fgxf^ {Z34 + Zi2{dUidU2 + dU3dU4)} |1, 1 >5 

/lo|0,2 >=/io{2;i2^34 + \{\Z12\^ - \z34\^){duidU2 + dU3dU4) 

^ (3.17) 

— Zi2Z34dUidU2dU3dU4^{dv\dV2 + dV3dV4)\i) > 
fllX^'\0,2>f'=fnx''{zi2Z34+li\zi2\'^ -\z34\'^)idUidU2 + dU3dU4) 

^ (3.18) 

— Z\2Z34dU\dU2dU3dU4\\l >5^ . 

Our choice of normalization in constructing these forms is described in Appendix C. We 
take the ground state to be the sum of these eleven forms. 
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3.3. Dynamical constraints 

What remains is to determine the consequences of the eight equations, 

Qa* = 0. (3.19) 

First note that each term in Qa can be assigned a parity (±, ±) according to whether it 
changes the parity of the wavefunction in {x, q) respectively. For example, the term 

i^^s^s^^)ax^'qj (3.20) 

has parity (— , — ) since it is odd in x and odd in q. Likewise, each term in \1/ has a definite 
parity. We can therefore isolate all terms in with a definite parity. It is also sufficient 
to restrict to the case a = 1 because our wavefunction ^ is Sp{2) invariant, but an Sp{2) 
transformation rotates us from one choice of charge to another. 

We can then ask: what combinations give terms with parity (+, — )? A quick check of 
Qa acting on the possible forms composing gives the following equation, 

{s'^)aPj { /i|0, > +/io|0, 2 > +/2IO, 4 > +/3I2, 2 >i +hx^x''\2, 2 } + 
(7Vs2^)„xV{/4a^l2,2 >g +/na:^|0,2 (3.21) 

+ {r>^)aP'' {fsX^ll, 3 >P +f9XP\l, 1 >^} + {DX)a {Mh 3 > +Ml, 1 >} = 0. 

The terms giving (— , — ) satisfy: 

(7VsV)axV{ /i|0,0 > +/io|0,2 > +/2|0,4 > +/3|2,2 >i + 

75^^2,2 } + {s^i;)aPj{hxP\2,2 >g +/nx^|0, 2 >^} (3.22) 
+ iDX)a {fsx'^ll, 3 >P +fQxP\l, 1 >f} + irX)aP'' {/ell, 3 > +Ml, 1 >} = 0. 
From (— , +), we find: 

(7''A)ap^ { /i|0, > +/io|0, 2 > +/2IO, 4 > +/3I2, 2 >i +hx^x''\2, 2 } + 
{DX)a{hx''\2,2>P +f^^xP\0,2>P} + {s^^|;)aPj{f8xP\l,3 >" +hxP\l,l>P} (3.23) 
+ (7^5^^;), x^q^ {/ell, 3 > +M1, 1 >} = 0. 
The last equation follows from considering the (+, +) terms, 

(DA)4/i|0,0 > +/io|0,2 > +/2|0,4 > +/3|2,2 >i +/5a;^a;'^|2, 2 }+ 
{rX)aP^{hx''% 2 >g +/iix''|0, 2 >P } + 

3.24 

{Ys^s''i;)ax''q^{hx''\l,2>>P +/9x''|l, 1 } + 

(s^V)aP,{/6|l,3>+/7|l,l>} = 0. 
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In each equation, we set a = 1 as a first simpfification. After evaluating angular derivatives, 
we are free to rotate q using Sp{l)ii so that ^ ^ and q-* = for z > 1. In a similar way, 
we can consider the point ^ Q with x'^ = for > 1 after evaluating the x angular 
derivatives. With this choice of coordinates, y = \q^\ and r = \x^\. 

The first set of equations relate fi^fr^fQ., fio and /n. These follow from considering 
the (4, 1) forms and the (3,0) forms in (3.21), (3.22) and (3.23), (3.24) respectively: 

r^-y§|-^/r + 5/,-4/, + 2/,o = 0, 
^^'+V|A-#U2r/n=0, 



r^- + y^- + ryfr = 0, 
oy or 

Likewise, by considering the (0,3) forms in (3.21), (3.22) and the (1,4) forms in (3.23), 
(3.24), we find the following equations relating /2, /e, /s, /lo and /n: 

+ y^ + ^/e + 5/8 + 4/2 - 2/10 = 0, 



dr " dy 2 
^■^'+V<!/2 + ^K2r/ii=0 



(3.26) 



dr " 2 

df2 9/8 

r^- + ryfe = 0, 

or oy 

9f6 2 f .y^ f n 

^-ry/8 + y/2 = 0. 
Note that equations (3.26) are the same as (3.25) under the identification: 

/2^/l /6^/7 /8^-/9. (3.27) 

It is easy to check there are no non-vanishing (0, 1), (4,3) and (1,0), (3,4) forms in 
(3.21), (3.22) and (3.23), (3.24). This leaves equations coming from forms with degree 
(2,1), (2, 3) in (3.21) and (3.22), and (1,2), (3, 2) in (3.23) and (3.24). Let us start by 
considering the (2, 1) parts of (3.21) which give the following constraints, 

/lo + 77^^ 2/1 +4/9 + r /4 + - + -T—) + TT^ /ii = 0' 

2 oy y oy oy 2 



f .ydfio 2f . dfg 1 a/3 4 2^/5 

/10 + 77^ 2/1 -r /4 + r— + /9 - - — + -r — ) 

2 oy or y oy oy 

+ %h + '^r^hi = 0, 

7 — ^ r— 3/9 + — /7 = 0. 

b y ay y ay or 2 



(3.28) 
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The (2, 3) terms give the equations, 

f ^ydfio 9/8 2. 1.9/3^429/5. 

/lo+o^ 2/2 -r^ /s-r /4--h^ + -r — ) 

2 9?/ 9r y dy b dy ^g) 

+ y /6 + yrVii = 0, 

2r^dh 29/3^ ^/s^o, n 
b y oy y ay or I 

Note that equations (3.29) are the same as (3.28) under the identification (3.27). This is 
a nice check that the equations are correct. 

We next need the (2, 1) parts of (3.22) which give the relations, 



, ,^42^19/4, y9/n 



^/lo + /3 + -r'h + -^ + 5^ + /ii = 
2 5 y ay 2 ay 

f f ^ 2f 1 9/4 ^ 1 9/7 2/% 2/ 9/11 . r (o on^ 

2 5 y ay r or 2 2 ay 

Of 19/7 

From the (2, 3) components, we find the equations: 

f , ^ , 4 2 . . 1 9/4 2/ 9/11 
^/lo + /s + -r- /s + + + hi = 0, 

2 5 y ay 2 ay 

y% . 4 2. 19/4 19/6 y^f.ydfii 

^/lo -/3--r/5---^ + --K irfs + o -K— + hi = 0, (<5-^l) 

2 5 y oy r or 2 2 oy 

Of 2 2. 19/6 y% ^ 

2/3 - -r h - ^78 = 0. 

5 r or 2 

These equations are again consistent with (3.27). 

We now turn to the (1, 2) parts of (3.23) which imply that, 

9/10 ^19/3 4 9/5^28 1 2, ^2/% . 

+ + + ^/s + o2/ /4 + 2/9 + — /ii = 0, 



2r dr r dr b dr b 2 4 
2/29/10 19/3 49/5 28 1 2, ^ 2f 

n — a s — rh + T.y U + y h 

2r ar r or b or b 2 

+ y-^ + 2/8 - -r/11 = 0, 
dy 4 

29/3 2 9/5 14 2, ^ ^/s^of Of n 

— T'':^ r/s -y f6 + y-^ + 2/8 - 2/9 = 0. 

r or b or b ay 
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(3.32) 



The (3, 2) forms give the following set of equations: 

a/io , 1 a/3 , 4 a/,. 28 1 2 Of n 
7^'^ + + 7;y h- 2/8 + -r/ii = 



2r dr r dr 5 dr 5 2 4 

2/2a/io ia/3 4 9/5 28 ^ 1 2^ , 2^ a/< 



f/s + T^y /4 + y /? - y- 
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2r dr r dr 5 dr 5 2 dy 

2/9 - ^/ii = 0, 

29/3 2 a/5 14 2, 5/9 . 

^r- -f5 -y fr- 2/9 + 2/3 = 0. 

r ar o or 5 ay 

Again, (3.33) and (3.32) are identical under (3.27). 

The (1,2) parts of (3.24) give the following equations, 

-r/10 + 17/3 + -^—h + + 5/4 + 2/7 + — /ii + — — — = 0, 
4 2 5 or 2 2 or 

-J ho - y /3 - ^/5 + + 5/4 + + 2/6 

^22. 2/% y^T^ 5/11 n 
+ ry /8-y/n-^^ = 0, 

//a - ^/5 - - 2/6 + 2fr + r'y'fs + y'h^ = 0. 
5 dy 

Lastly, the (3,2) parts give the equations: 



y^ 2y^r^ 5/4 y^ y^r 9/ 



-/lo + ^/3 + ^^/5 + + 5/4 + 2/6 + ^/n + 
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4 ■ 2 ■ 5 ar 2 2 5r 

y\ y\ '2y'r\ dh 
2 2f y^ f y^r dfii 

-ry /9-y/n-^^ = 0, 

y'fs - ^/5 -y^ + 2/6 - 2/r - r2yV9 + y'fii = 0. 
5 dy 

Again, note that (3.34) and (3.35) are identical under the exchange (3.27). 



(3.33) 



(3.34) 



(3.35) 



4. The Structure of the Bound State Wavefunction 

4-1- Reducing the number of functions 

Initially, we have eleven independent functions obeying a set of coupled differential 
equations. To make progress, we need to reduce the number of functions in a systematic 
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fashion. We can begin to whittle down the number of independent functions in the following 
way: the difference of the third equations of (3.30) and (3.31) together with (3.25) imply- 
that, 

/i = /2. (4.1) 

Let us now turn to (3.25). By taking dy of the second equation and dr of the fourth 
equation, we find two equations for d^^fr. For these equations to be compatible, we 



require that: 



/io-/i + /9 + -^ = 0. (4.2) 
y oy 



This equation is actually not a new addition to our list of constraints. It follows from 
(3.30) and (3.25). The same analysis applied to (3.26) gives, 

/io-/i-/8 + -^ = 0. (4.3) 

y dy 

These equations together require that, 

/s = -/9- (4.4) 

The difference of the fourth equations in (3.25) and (3.26) then implies the equivalence, 

fe = /7. (4.5) 

In this way, we are reduced to eight functions {/i, /s, /4, /s, A, /g, /lo, /ii}. Half the 
equations we derived are now redundant since the symmetry (3.27) is an actual identity. 

We obtain an algebraic relation between the remaining functions by using the third 
equation in (3.30) together with (3.25): 

2 2 

/3 + y /i - y /5 + fii = 0. (4.6) 

Using the algebraic relation, we can eliminate one function. The remaining 7 functions of 
2 variables are constrained by 14 equations, which is the minimal number we could have 
expected. 
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4-2. A first reduction of the equations 

Let us summarize the equations that describe the bound state. We have taken some 
simpler hnear combinations of the previous equations: 



(1 
(2: 

(3: 

(4 

(5, 
(6 
(7 

(9 
(10 

(11 
(12 

(13 
(14 



— — + 5/9 - 4/1 + 2/10 = 0, 

or ay 2 

+ryf7 = 0, 

oy or 

y/9 + y/i = 0, 

2/10 + - 4/1 + 5/9 + y r /ii + r— + y = 0, 



y 



2r^|//4 - + 32//9 - ry^ + 2-^ + -r 



5r 



9y 5 dy 



= 0, 



2r2a^ _ 2a^_^a^ 
5 y dy y dy 



dr 



3/9 + y /7 = 0, 

y fio + -Trh + 2/11 + y^ h = 0, 

2 oy r or 

2y/3 + -r yh - — /g + 2- = 0' 

5 2 oy r or 

^y(/4 + /t) - r— + y-^ = 0, 

-r^'/s - ^^1/ /7 + 4r/9 + -ry /n + ry— + 2— + — 
5 2 oy or b or 

2dfs 2 a/5 14 2, n 

— -y 17- y^ — 4/9 = 0, 

r or o or o oy 



= 0, 



10/4 + 4/7 - rVh + ^/lo + + 2r^ = 0, 



9/7 



a/4 



dy 



yh + t^^yh + r^yfd + yfn - ^p- 

5 dy 



+ ry 



dr 
df 



11 



dr 



= 0. 



(4.7) 



In these equations, we can remove one function using (4.6). This is a comphcated set of 
coupled equations. To uncover the structure of the wavefunction, let us begin by simplifying 
as much as possible. 

After staring at these equations for sometime, a pair of equations (3) and (4) in (4.7) 
appear distinguished. At our special point, these equations involve only the (4, 0) and (3, 1) 
forms given explicitly in Appendix C. They do not involve any (2, 2) forms. These are the 
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analogues of the 'top forms' which played a crucial role in proving non-renormalization 
theorems M. We can eliminate /g from these two equations giving, 



[-di + ^dy + rV] fr = y'{l + \dy - rdr] /i, 



{A.i 



Note that if we set the source term 5'(/i) = 0, then the homogeneous solutions for fj are 
g±ry /2 rpj^g pj^g sign is not normalizable. However, both solutions suffer from a more 
serious problem. The ground state wavefunction must be a smooth function. It must have 
a convergent Taylor series about the origin. This implies that each fi must be function of 
and near the origin. The homogeneous solutions alone are therefore ruled out. 
What this teaches us is that /y is determined in terms of /i. It is convenient to make 
the following redefinition, 

/r = 2/77- 

Note that /y can have a term near the origin. The equation (4.8) now takes the form. 



-dl-hy + rW\f7 = S{h). (4.9) 



The left hand side of (4.9) is the Hamiltonian for a four-dimensional harmonic oscillator! 
Now there is a very pretty collapse. Equation (4) in (4.7) determines /g in terms of /i. 
Likewise (1) and (2) determine /lo and /n. Equation (10) determines /4, while (14) 
determines /a. The algebraic constraint (4.6) then fixes /s. Clearly, there many other 
ways to collapse the problem. The main point is that all the remaining functions are given 
in terms of /i. This leaves 8 equations which must determine /i. 

4-3. Solving for f J 

We can now express fi in terms of /i in the following way: 

/r = /r°(r2)e-^^'/^ + {-^J - hy + r^^^ ' S{h). (4.10) 

The y = component is determined by requiring that /y be smooth, as we discussed 
previously. Smoothness of /y also requires that be a function of r^. In turn, we can 
expand /i as follows: 

oo 

/i = E/i''WI^>- (4.11) 

n=0 
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The \n > are radial harmonic oscillator eigenstates which obey, 

- ^dy + rVj 1^ >= En\n >, 

where En = 4(n+ l)r. The construction and properties of these eigenstates are described 
in Appendix D. These eigenstates have the nice feature that acting on \n >, the operators 
dr, ydy and involve only \n — 1 >,|n>,|n+l >. Using the relations from Appendix D, 
we see that the source term has a beautifully simple form: 

Sih) = J2ifi-r9rf^)\^>- (4-12) 

n 

It is now easy to solve for /y in terms of \n >, 

oo ^ 

fr = f^{r')\0 >+y'J2-^ - "-^rfi) \n > . (4.13) 

We have left in (4.13) for later convenience. By considering the coefficient of y^ in (4.13) 
and imposing smoothness, we obtain the following relation: 

oo 

-i/7 + E:^(i-^^-)/i" = o- (4.14) 

n=0 

Since is non-singular as r— >0, we obtain the sum rule 

oo ^ 
n=0 ^ 

and the relation: 

9 °° 1 

n=0 ^ ' 

Note that these sum rules may be largely formal since we do not know whether the sums 
are absolutely convergent. 

4-4- Equations for the physics near the flat directions 

Instead of proceeding to reduce the number of functions, let us take a different tack. 
The most interesting physics in this problem occurs in a neighbourhood of the flat direc- 
tions. So let us consider a Taylor expansion about the flat direction y = 0. This approach 
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turns out to be more useful than reducing the number of functions, which increases the 
complexity of the resulting equations. We expand each /j in the following way, 

/.-t°(r')+^?(rV+tt(rV + ---- (4.17) 
The algebraic constraint (4.6) together with the equations of (4.7) give the following set 
of relations on the t^: 

(1) tl-jtl + t^,, = o, 

(2) rdr4 + 5^9 - 4t? + 2t% = 0, 

(3) drt*^ + 2rt?i = 0, 

56 8 

(4) —r4 + 4r4 + 2dr4 + -^r^dr^ = 0, (4.18) 

(5) ldr4 - Irdrtl - y - 4tg = 0, 

(6) ml + 4^ + 2rdrtl = 0, 

(7) tl + 24 + \{t'i+t%] = Q. 

This gives 7 equations for 8 unknown functions. As we might have expected, this is not 
sufficient to determine the flat direction physics without input from higher y terms. We 
can similarly derive equations involving only tf and t° which are given in Appendix E. 

Note the critical observation that we can solve for all t| in terms of the t° using 8 of 
the equations of E.l. There are no new independent functions at order in the Taylor 
expansion. We can express the in terms of the t° and at most their second derivatives 
using the first 7 equations and equation (10) of E.l, 

(1) t'7--\rhl 

(2) ti = -^Kt? + rt?}, 

(3) tl^^\{'itl+rdrtl-2tl], 

(4) tl = -^^{Ttl-drtl], 

(5) tl=^{-2rHl-1^4-Zrdrtl], 

(6) tl = -^{l0rt° + 8r3t0-5a,t?}, 

(7) = 4 {~^rH\i + -dA + rdrA + dl4 

lb y r 

(8) tl = --^ {8rt? + 2rhl, + m'l + r'dr^ + rd^,^} . 

lor 
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(4.19) 



What we need to check is whether any of the remaining 7 equations of E.l are new 
relations. After some algebra that we will spare the reader, it turns out that all the 
remaining equations of E.l are consequences of (4.18). This is not too surprising. The 
interactions involve order terms so we should not be able to completely determine the 
physics on the flat directions without expanding to higher order in y. So let us expand 
to order giving equations relating the t|,tf and t^. These equations are again given 
explicitly in Appendix E. 

The equations of E.2 are very similar to E.l, except for mixing with some terms 
through the and interactions in (4.7). It is easy to see that we can again solve for 
all the tf in terms of the t^. It should be clear that all the coefficients of the higher y^" 
terms in the Taylor expansion are determined in terms of the t^. The problem is then 
to determine the t^, and we need one more relation in addition to those of (4.18). Again 
using the first 7 equations and equation (10) of E.2, we can solve for the tf: 

If 1 

— <8tl + 2rdrtl - 8tl + -drtl 
lo { r 

^ {rt? - 4rHl + 2rtl + drt^ + Ar^drtl} , 

^ {rt? - 4rhl - 38rtl - 4c»,t? - Gr^drtl} , (4-20) 



(1) 


- - 


(2) 


^9 — 


(3) 


''11 — 


(4) 


^5 — 


(5) 


- 


(6) 


- ■ 
^4. — 


(7) 


^10 — 


(8) 


ti = - 



200 I " " ' r 

1 



{ - 5rH^ - lOrHl + 2rHl + 45rt^ + 2Qdrtl 



200r 

Note that the expression for agrees with the expression coming from the earlier relation 
(4.8) that we derived between /i and /y. Once again we are left with 7 additional equations. 
It turns out that these additional equations again give no new relations. After checking 
higher order Taylor coefficients, we find that there are no further relations. It appears 
that any choice of t° satisfying (4.18) give a zero energy solution. However, most of these 
solutions are not normalizable. 
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It is straightforward to derive the general recursion relation for tf in terms of lower 
Taylor coefficients: 

(1) t:} = -l\rHr'+ltr'], 

(2) t- = -^{a,t--2 + rt^-2}, 

nr ^ ' 

(3) ^ri = ^{^^9-'-^r'- J^r^? 
1 1 1 ^ .„_2 ,n-2 



(4) t^ = -|-a,tr -^4 

(5) n-V^n-nx-\tr^ 



(4.21) 



(6) ^l-Yn{ J^^'""' + - 2'^" V'^" 

(7) t^^o = - 2n{^^ 6) ^^"^ + + ^^''^^ + + + 2nra,t^} , 

(8) t- = -l{t--2+nt^o+r'iir'}- 

Let us close this discussion of the Taylor expansion by pointing out that the entire Taylor 
series depends only on 3 of the t\. To see this, let us express (4.18) in a form which will 
be more convenient for later manipulation, 

(1) ^^ = ^(^^+^?i)> 

(2) t\-^ = —-^drt^, 

(3) t'^ = -^(5tUrdrtl), 

^9 = 6r^''^^' (4.22) 

(5) tiQ = 2ti ~ ~ 2^rtg, 

(6) 20dr{rh'^s) + ^ {729^ + ^^rd^ + ^r^d^}4 = 0, 

(7) drtl + ^^{-64r^ + 72dr - IGr^dr - 72rd^ 
+ 12r^d^ + 12r^d^ + r% } = 0. 

The first 5 equations of (4.22) are definitions of various The last 2 are relations on the 
3 independent functions t\,t% and ^4. 
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5. An Asymptotic Expansion of the Wavefunction 

5.1. Matching Taylor and oscillator expansions 

It is hard to see how to implement the normalizability condition in a Taylor series. 
So although we have found the zero energy solution in terms of the t^, we need a practical 
procedure to construct the t^. We will determine the t^ under the assumption that the 
bound state wavefunction admit an asymptotic expansion in powers of 1/r. The primary 
motivation for studying the asymptotic expansion is that the asymptotic form should 
be interpretable in terms of supergravity plus higher derivative corrections. What we 
do not know is whether the asymptotic expansion converges to the actual bound state 
wavefunction. This issue is closely related to the following two questions: if we sum the 
effects of all higher derivative terms in the effective action for this gauge theory, does the 
result give non-singular physics as r^O? If we sum the effects of all higher derivative 
terms beyond supergravity on the spacetime solution for a 5-brane in M theory, do we find 
a smooth convergent solution near the 5-brane? 

These questions are intrinsically tied to the problem of gauge-fixing the Taylor series 
solution in a way that results in a normalizable solution. We want to construct the t^ 
in a useful systematic expansion. The dominant terms in an asymptotic expansion are 
those that decay polynomially in 1/r. We point out that for large r, an approximate 
asymptotic bound state can be constructed in a 1/r expansion using a method described 
T^. An analytic expansion in 1/r near infinity is essentially a perturbative expansion. 



m 



although not in a priori but in g'^^^ . A similar technique was used in |1T5| , |T6i|T7| to 



further explore the long distance dynamics and the asymptotic structure of the bound state 
wavefunction for 2 DO-branes. Unfortunately, the effective long distance Hamiltonian has 
only been constructed to order 1/r^, which is the required order for an index computation 
0,1181. With our knowledge of the Taylor series solution (4.21) for this problem, we can 
do significantly better than those approximate constructions. 

Can there be non-perturbative terms? These are terms which are not visible in a 

2 

1/r expansion, like , but which become important as r becomes small. We do not 
actually know whether there are any such terms, and there are no candidates like instanton 
configurations that could generate these terms in the abelian gauge theory. This leads the 
first author to suspect that the analytic expansion in 1/r might well be exact. Nevertheless, 
we cannot prove that non-perturbative terms are not present. To really rule out such terms, 
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we need a technique for finding the bound state solution which is inherently more global 
than an asymptotic expansion. We shall discuss a more global approach in section 6. 

What is hard to see in the Taylor expansion of section 4 is normalizability in the 
y-direction. This is much easier to see in the oscillator expansion so we will match the 
two expansions to get control over the question of normalizability. When expanded in a 
harmonic oscillator basis with frequency proportional to r, 

oo 
n=0 

the fl^ must decay as r—>-oo. Let us expand each fl^ in powers of 1/r. This implies that 
the have an expansion in 1/r. We can then reorganize the Taylor series for fi in the 
following way: 

p k 

The 6^'* are just some numbers which determine the collection of harmonics contributing 
to a given power 1/r^. Note that the oscillator eigenstates depend only on the combination 
ry^ so that '^j^ ^'*|^ > is a power series in ry^. 

For example, suppose that the harmonic |m > is the only harmonic with a non-zero 
coefficient in the sum ^'*|^ > ^or the 1/r^ term of fi. We stress that in general there 
can be many harmonics contributing to a given term, but for simplicity, let us assume 
there is just one. The Taylor series for fi must then contain the terms: 

„. 1 



= {^1^^ - r } + {^^"^ - r^^"^^ + i"'} + •••)+••• ' (5-3) 

= 6^^ (^1 - ^(1 + mW + ^ (3 + 4m + 2m')r^y'' + . . + . . . . 

There are specific relations between the Taylor coefficients in (5.3). We want to impose 
relations of this kind on the ti to satisfy normalizability in the y-direction for each choice 
of p in the 1/r expansion. 
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5.2. The structure of the solution 

Our Taylor series solution is completely determined by the three functions ^1,^3 and 
^4. These three functions must obey equations (6) and (7) of (4.22). To ensure that t1 is 
decaying, we see from (6) that t'^ must take the form, 

*§ = ^ + ..., (5.4) 

where omitted terms decay more rapidly. For a given choice of p in the 1/r expansion, we 
wish to extract the terms in each ^3 which contribute to ^fc'^l^ >• Again, this is just 
the statement that we can organize the Taylor series for fs so that, 

/3 = ^E^^''l^ >+•••' (5-5) 

k 

3 3 

for some coefficients 6^' . This is formally true for any Taylor series. What is generally not 

3 3 

true is that generic 6^' give a wavefunction normalizable in the y direction. Heuristically, 
the norm of should be dominated by the leading term in the 1/r expansion. If we 
compute the norm of /3 under this assumption, we see that: 

where we have taken the normalization of |A; > given in Appendix D. 

By looking at (4.19), (4.20), (4.21) and (4.22), we can see what terms in the Taylor 
series for have the right structure to determine the 6^'^ coefficients of (5.5). To quickly 
answer this question, let us list some order of magnitude relations which follow from (4.22) 
for the perturbative expansion. All relations are given in terms of our 3 independent 
functions t^,t^ and t^. 



(1) 


tir.o{tiy + tyr% 


(2) 




- 0{t\lr\ 


(3) 


h ' 




(4) 


h ' 


- 0{tl/r% 


(5) 


''10 


~0(t° + t°/r^). 
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There are similar order of magnitude relations for tf from (4.19), 

(1) t'7-o{4), 

(2) tlr.Oityr' + t'i), 

(3) t?,~0(t°/r^ + t°), 

(4) tl^O{tU4/r'), 

(5) 4^0{rX + 4/r% ^^'^^ 

(6) tlr.O{4 + 4/r'), 

(7) tl,r.0i4 + tyr'), 

(8) 4^0(4 + 4^). 

It is easy to continue and list order of magnitude relations for tf from (4.20), and higher 
tf using (4.21). These relations are useful for easily determining which terms in tf are 
relevant for determining the oscillator coefficients 6^'*. 

Returning to our specific case of /s, we see that since t^ ~ ^/r^, the only terms in t^ 
relevant for computing 6^'^ are those proportional to 1/r^. Looking at t^ from (5.7), we 
see that the only way to have a term contributing to the ^ b^'^\k > is if. 



t2=^ + .... (8.8) 



It could be the case that c[ = 0, which means that the 6^'^ sum up in such a way that the 
y^/r^ term in (5.5) vanishes. Actually if c[ = 0, the situation is much worse: we can see 
after some work that the 6^'^ have to be chosen so that the (ry^)^"*"^" terms in the Taylor 
expansion of. 



k 

,2\2n 



vanish for all n. The {ry ) terms are all proportional to ci and give the constraint. 

However, the left hand side of (5.9) is the expansion of the lowest oscillator |0 > with the 
^^y2Y+2n ^erms missing. Since |0 >= e"''^^/^, we immediately see that the left hand side 
of (5.9) sums to the expression: 
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This is clearly not normalizable so is necessarily non-zero. We will use this kind of 
argument repeatedly to determine the asymptotic solution. 
With 7^ 0, we get a modified constraint: 

Cl (l + ^.ry'f + ...)- c; ("^(ry^) + ^^{ry'f + ...") = |fc > • (5.10) 

It is easy to see that the terms proportional to sum up to give,i 

10 I 

In this case, it is completely clear that we only get a normalizable solution if we pick, 

c'l = 5ci. (5.11) 
More generally, we will encounter the situation where we have a Taylor series of the form, 

aY^d^ {ry'f- + a'Y,d'r. iry'r''+\ (5.12) 

n n 

where we know both and d'^, and either a or a'. Let us assume we know a and we wish 
to determine the values of a' for which the Taylor series is normalizable. If the Taylor 
series can be fitted by a finite number of oscillators for some a' then that choice of a' 
is unique. To see this, suppose there were two distinct choices of a' . We could take the 
difference between the two series to obtain a Taylor series proportional to, 

n 

However, this series cannot be generated by any sum over a finite number of oscillators 
since all the (ry^)^" terms must vanish. Therefore the choice of a' is unique. 

From equation (6) of (4.22), we see that the 0(l/r^) term in induces an 0(l/r^) 
correction to t^. This begins to suggest a perturbation expansion in 1/r^, and indeed that 
seems to be the case. From t\ ~ 0(l/r^), we see that we need a 0(l/r^) correction to t\ 
and so on. It is easy to add t\ into the story: from the expression for t\ ~ 0(l/r^), we 



To see this, let us introduce new variables {r,u = ry^ /2). In terms of these variables and 
simple redefinitions of the /i, the equations of (4.7) can be put in a triangular form with respect 
to the grading induced by the 1/r expansion. In doing this, we treat u as independent of r. The 
resulting triangular system is exactly solvable in the 1/r expansion. 
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learn that ~ at leading order. From (7) of (4.22), we see that t° mixes with the 
first correction to and so on. What we have uncovered is the minimal required form 
for the solution given ci 7^ 0. The choice of ci then determines the normalization of the 
wavefunction. The solution is given by the expansion: 

^ ' 

^4 = ^ + ^ + ^ + ---, (5.13) 

^" Jl r" 

It is interesting that the have an expansion in powers of Restoring the coupling 

constant, we see that this is an expansion in jr^ . This is the natural expansion parameter 
in the gauge theory. From the perspective of the effective action on the Coulomb branch, 
what we are doing is summing the effects of the metric and all higher derivative corrections 
on the vacuum state of the gauge theory. 

The expansion (5.13) represents the minimal required terms in the solution. It is 
natural to ask whether other powers of 1/r are possible. Equation (6) of (4.22) requires 
that, 

except for the special case where ^3 ~ 0(l/r^). Likewise, equation (7) of (4.22) requires 
that, 

t\ ~ 0{rH\ + lt°), 

again except for the case ^4 ~ 0{l/r'^). Let us be concrete: suppose ^3 has a term of 
order 0{l/r^). This requires ~ 0(l/r^) and therefore ~ 0(1), which is not possible. 
Suppose has a term of order 0(l/r^), which implies that tjj ~ 0{l/r^). In turn, we see 
that t| ~ 0(l/r) which requires that ^3 = ajr^ . This is impossible unless a = 0. In that 
case, the Taylor expansion of jj, has terms of the schematic form: 

/3-0(^) + 0(r/) + .... 
r 

As in our earlier discussion, it is not hard to check that this sum is not normalizable. This 
kind of argument extends to higher powers: suppose t\ has a term of order 0{l/r'^), which 
implies that t% ~ 0(l/r^). From ~ 0(l/r^), we require a 0(l/r^) correction to which 
brings us back to an earlier case. This kind of reasoning suggests that only those terms 
generated by the original 1 /r^ term in appear in the solution. This also agrees with our 
gauge theory intuition. We shall therefore restrict our attention to the terms of (5.13). 
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5. 3. Supergravity and the leading terms of the solution 

So far, we have found that the leading order term of /a is proportional purely to the 
oscillator ground state |0 >. We also know all the leading powers for the t^. A glance at 
Appendix C tells us that the dominant terms for large r in the bound state \1/ are /a, 
and /s- Note that on the flat directions where y = 0, the only non- vanishing forms in \E' 
are those with coefficients /a, and /s. It is not hard to check that both and /s are 
also proportional to |0 > at leading order. Recalling that r = |a;^|, we have learnt that 
these leading terms are given by, 

/3 = ^|0>+---, 

xV4 = ^|0 >+..., (5.14) 



,l^2 . 



(xT/5 = ^^#^|o > + .... 



Although each term transforms very differently under Spin{5) in either the 1, 5 or 14, 
each decays at the same rate at leading order. 

We can compare what we have learned about the exact Taylor series solution with what 
we might expect from an effective Hamiltonian construction like the one given in . To 
construct the approximate bound state, we note that the wavefunction is sharply localized 
near the flat directions for large r. The potential term (2.20) is dominated by the r^y^ term. 
We conclude that for large r, the wavefunction can be expanded in a harmonic oscillator 
basis, with higher oscillator modes suppressed by powers of 1/r. The approximate bound 
state wavefunction takes the form of a product: there is a wavefunction of the light vector 
multiplet degrees of freedom multiplied by the ground state for the massive hypermultiplet 
degrees of freedom. We can now construct the effective Hamiltonian governing the long 
distance physics in a 1/r expansion. The hypermultiplet ground state is easily determined 
from (2.16) to have the form,0 

re-''y^/^dvidv2. (5.15) 

In the most primitive approximation where the effective interactions are cancelled to order 
1/r, the effective Hamiltonian is {l/2)p^p^ and acts on (5.15) multiplied by some function 



^ The natural decomposition of and qi into light and heavy variables at large r makes the 
construction of the approximate ground state much simpler than the non-abelian DO-DO theory 



considered in |14]. 
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of X and du. In fact, to determine the leading decay, we should go to order Even 
then there will be a degenerate set of approximate ground states. 

However, this rough approximation is good enough for the purpose of comparison 
with the structure following from the Taylor series solution. Using the forms given in 
Appendix C, we see that the coefficients in (5.14) precisely conspire to give agreement 
with the approximate construction. The leading terms in ^' sum to give, 

^ = fs\2, 2 >i +Ux^\2, 2 >l +h{x^)^\2, 2>\\+..., 

lOci 2/ (5.16) 
~ — —e~^^ ''^{duidu2 — dusdu4){dvidv2) + . . . , 

where we take x^ positive. There are a number of strange features of this asymptotic 
solution. First note that the vacuum for the massive fermions dvidv2 is not a single repre- 
sentation of Spin{5). Let us contrast this with the case of 2 DO-branes where the vacuum 
for the massive fermions transforms in a single representation, the 44, of Spin{9) |lT5|,p!7 



The appearance of three different representations, the 1,5 and 14, in the asymptotic so- 
lution really cries out for an interpretation both in terms of the DLCQ M theory 5-brane, 
and in terms of the supergravity solution for the D0-D4 bound state. 

What this seems to suggest is that the massive degrees of freedom never really decouple 
at large r. Through the Spin{5) flavor symmetry, the vacuum always knows about the 
existence of massive degrees of freedom. This issue is intimately tied to uniqueness of 
the bound state: the statement that the bound state is unique involves knowledge about 
both long and short distance physics. Uniqueness however forces invariance under the 
full Spin{5) acting on both light and heavy degrees of freedom. That dvidv2 is not an 
irreducible representation then requires particular combinations of Spin{5) representations 
for the light degrees of freedom. A deeper understanding of this issue is certainly in order. 

5.4- Beyond supergravity 

Let us return to our general solution (5.13). We determined c'l in terms of ci in (5.11). 
How do we determine ? The straightforward way to determine is not particularly 
elegant. We can take the Taylor series for fi and impose normalizability in the y-direction. 



r^ 



J' 



(5.17) 
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Note the general structure: half the terms in the Taylor series are known. The other half 
depend on the unknown constant that we wish to determine. It is not hard to see that 
the (ry^)-'^"'"^"' terms come from expanding the lowest oscillator |0 > state. The situation 
is then essentially the same as in (5.10), and we find that: 

c'l = —c,. (5.18) 

From equation (6) of (4.22), we find that C2 is determined by c[: 

C2 = ci. (5.19) 

Equation (7) of (4.22) relates C2 to c'/, 

/ 25 , , 

c'2 = -yci. (5.20) 

We are again back to the question of studying /i to determine C2. The relevant terms in 
/i begin at order 1/r^, 



(5.21) 



Apart from more complicated coefficients, the general pattern is the same. Half the terms 
in the Taylor series are proportional to ci while the other half depend on the unknown 
constant C2. We can try to match (5.21) with an oscillator expansion involving a finite 
number of oscillators. Indeed, for the choice 

1 7^ 

4 = -^ci, (5.22) 

we can fit (5.21) by the first two oscillators: 

225 ci /, 2, \ 
/i = -^Z^ |0>-n|l> +•••• 



8 r8 V 9 

As we might expect, higher excited oscillator states appear in the solution as we study 
more rapidly decreasing terms. This also suggests that only a finite number of oscillators 
will appear at any given order in the 1/r expansion. 
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Let us iterate the argument one more time. From (6) and (7) of (4.22), we immediately 
obtain the following relations: 



C3 = yci, 



35 
-( 

195 



4 = — ^ci. 



The relevant terms in fi take the form, 



1 / 2825C1 2.^jc'i 595ci\ 

425ci ^ _ 2\3 I /n/_,2\4 



and can be fit by the first three oscillator modes, 

^_525_ci_ 
~ 32 

for the choice. 



208,^ 


115, 


|2>-2j|l>- 


' 3 1' 







(5.23) 



(5.24) 



^^' = -^ci. (5.25) 



6. A Remarkable Reduction 

6.1. Prolongation 

Studying the bound state solution in an asymptotic series is analogous to studying the 
M5-brane in a derivative expansion around the supergravity solution. Ideally, we want a 
more powerful technique to solve the vacuum equations. The aim of this final section is to 
present a more global approach to solving the vacuum equations. Hopefully, this approach 
is closer to the method an M theorist might use to study the M5-brane. 

We shall present a surprising reduction of the long list of equations (4.7) to a single 
scalar elliptic equation. The equation takes the form, 

(^A + S- V + Vr) w = 0, (6.1) 

where A = + dy, the vector field B with components {Br, By) and the potential W 
are rational functions of r and y. The function u = u{r,y) is a particular combination 
of the fi. In order to explain how we reduce our first order system to a single second 
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order equation, let us first consider the inverse process: prolongation. To understand the 
procedure, let us begin with an illustrative example: we take an equation of the form. 



+ Fyy + UF^O. (6.2) 

We can define new functions p = and q — Fy with which we can 'prolong' our scalar 
second order equation into a system of first order equations: 

Px + Qy + UF^O, 

Fa, = p, (6.3) 
Fy = q. 

We also have a compatibility relation, 

F — F 

^ xy — ^ yxi 

which implies an additional fourth equation: 

Py-qx^ 0. (6.4) 

We can express these relations as a differential system, 

dF = pdx + qdy, 

dp = adx + hdy, (6-5) 
dq = bdx - (a + UF)dy. 

We could repeat the procedure and prolong again by adding the equations, 

da — cdx + edy, 

(6.6) 

db = edx + {{UF)a; - c)dy. 

Each time we prolong a system of equations like this, we add two new unknown functions. 
These functions are the unknown derivatives of the functions comprising the previous 
system of equations. 
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6.2. Deprolongation 

Now we would like to 'deprolong' our system of equations (4.7).li We shall see that the 
system (4.7) containing 7 independent functions can be obtained by prolonging a single 
second order equation three times. In order to deprolong (4.7), we write the equations as 
an exterior system as before. We choose 5 functions Fi which are linear combinations of 
the initial seven fi so that dFi is expressible algebraically in terms of the original seven. 
We keep the 5 equations defining the dFi and discard the remaining two equations. For 
example, for this first step, we can make the following choice: 

Fi = fi + /lo, 
i^2 = r2/5-2/'/i, 

So far the equations remain first order. We then iterate this procedure until we arrive at 
the first prolongation of the scalar equation. We then make the obvious substitution to 
transform the first prolongation into a second order scalar equation. 

Let us summarize the results of the deprolongation. We obtain the following equation: 



B B W 

Urr + Uyy + ^ Ur -\ ^ Uy + — U = 0. (6.7) 

rF yF F 



We define u in terms of the variables s = and t = and the functions 

hi f4:, /s, /?, fioi fu- 

u = {2s- t) h - {^ts^ - \t^s - V - 1^ /lo + {^ts^ - \t^s - ^-t^ + 1^ h 

It would be interesting to find a geometric interpretation for (6.7) - perhaps in terms of 
some bundle over either IR^ or IR^. Let us list the rational functions which appear in (6.7). 
For B^, we find the expression: 

Br = - 224t^s^ + 576^^5^ + 321^8"^ - 12t^s + 4f - 12t'^s 

(6.9) 

+ U4ts^ - 24t^ - 448sH - 48s^ + 36, 



We are especially grateful to Robert Bryant for suggesting and explaining this reduction to 

us. 
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while By and F are given by, 

By =212t^s^ - 192t^s^ - m^s^ + 264ts2 - 216^^5 + 68t^s + 42^^ - ist^ + 9 - 64s^t, 

F =9 + Ut^s^ - 8t^s^ - 6t^ + 2At^s - 112^^5^ + - 2Ats^ + MsH - At^s. 

_ (6.10) 
Lastly, the function W which determines the potential W takes the form: 

W = - 105^2 + SQsH'^ + 38t^ + 144st - m^s'^ + 2As'^ - lAAt'^s^ - 32t^s^ 

(6.11) 

+ 24t^s^ - 256s^t^ + 224s^t + 384s^t^ - - 76st^. 

Equation (6.7) is remarkably simple by comparision with (4.7). That the equations reduce 
this way opens up the possibility of answering a host of otherwise intractable questions. 
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Appendix A. Quaternions and Symplectic Groups 



We will summarize some useful relations between quaternions and symplectic groups. 
Let us label a basis for our quaternions by {1, /, J, K} where, 

= J^ = K^ = -1, UK = -1. 

A quaternion q can then be expanded in components 

q = q^ + Iq'^ + Jq^ + Kq\ 

The conjugate quaternion q has an expansion 

q = q^ -Iq^ - Jq^ - Kq"^. 

The symmetry group Sp{l)ji ~ SU{2)ji is the group of unit quaternions. Let A be a field 
transforming in the 2 of Sp{1)r. If we view Sp{1)r acting on A as right multiplication by 
a unit quaternion g then, 

A^A^. 

In this formalism, A is valued in the quaternions. Equivalently, we can expand A in 
components and express the action of g in the following way, 

where Qab implements right multiplication by the unit quaternion g. For example, right 
multiplication by / on g gives 

q^ql 
^q^I-q^ -q^K + q^J, 

which can be realized by the matrix 



r-R 



/O -1 0\ 

10 

1 

Vo -1 o/ 



(A.l) 



acting on q in the usual way qa ~^ ^ablb- The matrices and realize right multipli- 
cation by J, K while 1^ is the identity matrix: 



J 



R 








-1 


° \ 











-1 


1 











\0 


1 





/ 



K 



R 



1^0 -1^ 

10 
0-100 

Vi y 



(A.2) 
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We define operators in terms of {l^, J^, K^^ 




In a similar way, the group Sp{2) ~ Spin{5) is the group of quaternion-valued 2x2 
matrices with unit determinant. We will view Sp{2) as acting by left multiplication on a 
field in the defining representation. So an element U e Sp{2) acts in the following way: 



Equivalently, in terms of components 

Lastly, we can give an explicit form for the gamma matrices (2.3) in terms of quaternions: 

I fl 0\ 2 1\ 3 / / 











J ■ 





^ U -1 j ' ^ V 1 y ' ' \ -/ 



In turn, {/, J, K} can be expressed in terms of the Pauli matrices cr* 



loy v^o;' vo-i 

as 4 X 4 real anti-symmetric matrices: 



, _ / -ia^ ^ _ 
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Appendix B. Forms and Representations of Sp{2) 



Using the complexification (3.7), we obtain a set of Hermitian 4x4 matrices from 
those given in Appendix A: 



7 



1 

-1 



7 



f = 



~i(7^ 
3 



1 

1 











f = 






^ 




7^=( 


I — icT^ 






-la 




We also need the symplectic metric or charge conjugation matrix, 



C 



/ 1 0\ 

-10 

1 

V -1 o/ 



which implements complex conjugation: 

C7C = -7* 



(B.l) 



Our forms dua and dva transform in the 4 of Sp{2). The representation 4 decom- 
poses into 5® 1. As an example, we can explicitly construct the 1 from (1, 1) forms in the 
following way, 

duCdv, (B.2) 



while the 5 is given by: 



du^^C dv. 



(B.3) 



It is not hard to check that these combinations transform correctly. Lastly, we need to 
consider 5(8)5 = 1© 10 ©14 since the 1 and the 14 appear in (3.14). The 1 is given by 
the form, 

J^d'^l^Cdudv^Cdv, (B.4) 

while the 14 has components: 

du 7^'^ C du dv 7^) Cdv --6^"'J2 ^ ^ dv. (B.5) 
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Appendix C. A Word on Normalizations 

To fix the choice of normalizations, we Ust expHcitly the forms appearing in (3.12) 
through (3.18) at the special point where 7^ with = for n > 1 and qi 7^ with 
Qj = for J > 1. All forms act on the canonical vacuum |0 > which is omitted, 

/i|0, >= fi{qifduidu2du2,duA, (C.l) 

f2\Q,A>=f2{qifdvidv2dv^dv^, (C.2) 

/3I2, 2 >-i_=f:i{{duidu2 — duzdu^{dvidv2 — dv^dv^) — 
2{duidu4dv2dv3 + du2du3dvidv4)+ 

2{duidusdv2dv4^ + du2du4dvidvs)} , 
f4X^\2, 2 >\=f4X^{duidu2 — duzdu4){dvidv2 + dv^^dv^), 



f^{x^Y\2, 2 >^^=f^{x^)'^{^{duidu2 - dusdu4){dvidv2 - dv3dv4)+ 
2 

-{duidu4dv2dv2, + du2du3dvidv4)— 
o 

2 

- {duidu3dv2dv4 + du2du4dvidvs)} , 
5 

/ejljB >=fQqi{dvidv2 + dvzdv4){du\dv2 — du2dvi + du^dv4 — du4dv3), 
/rjl, 1 >=frqi{duidu2 + duzdu4){duidv2 — du2dvi + dusdv4 — du4dv^), 

fsx^\l, 3 >^ = fsx^ qi{dvidv2 + dv3dv4){duidv2 — du2dvi — dusdv4 + du4dv^), 
f^x^li^ 1 >^=fgx^qi{duidu2 + duzdu4){du\dv2 — du2dv\ — duzdv4 + du4dvz)i 

1 



(C.3) 



(C.4) 



(C.5) 



/io|0, 2 >=/io ^{qif'{duidu2 + du^du4){dvidv2 + dv^dv4), (C.6) 
/iia;^|0,2 >^—fiix^ -^{qiY {duidu2 + du^du4){dvidv2 — dvsdv4). (C-7) 
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Appendix D. The Four-Dimensional Radial Harmonic Oscillator 

We want to construct eigenstates for the radial four-dimensional simple harmonic 
oscillator which satisfy, 



-dy - ^dy + r^y^ \ \n >= En\n >, (D.l) 



where = 4(n + l)r. The easiest is the ground state: 

|o >= e-^y"/^. 

A general eigenstate takes the form, 

\n >= (l + ai")y2 ^ _ _ _ ^ a^^^k^"") e"'^^'/'- (D.2) 
It is not hard to check using the nice relation, 

(n) 

that each coefficient am is determined by the recursion relation. 



m{m + 1) 



where Oq'^ = 1. Note that these eigenstates are not normalized, but they are orthogonal 
when integrated with the measure y^dy. The norm of these eigenstates is given by the 
formula, 

<:n\n>= ir—^-—. (D.5) 
' r2 (2 + 2n) ^ ^ 

We will need to evaluate various operators acting on \n >. The nicest are the three 
operators y'^,ydy,dr. These operators raise and lower by at most one unit: 

ry^\n > = -n\n - 1 > +2{n + l)|n > -(n + 2)|n + 1 >, 
ydy\n > = -n\n - 1 > -2\n > +{n + 2)\n + 1 >, (D.6) 
2rdr\n > = -n\n - 1 > -2\n > +{n + 2)|n + 1 > . 

Note that ydy is equivalent to 2rdr when acting on \n >. 
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Appendix E. Equations from a Taylor Expansion 

We list explicitly the equations that follow from a Taylor expansion of the /j in the 
y-direction and which only involve and t'^, 

(1) l4-^-^tl + ti, + tl = o, 

(2) rdrtl - 6tl - ^t^ + 5tl + 2tfo = 0, 

(3) a,t? + 2rt?i + r|t?-^t°| = 0, 

(4) 2rtl + + rt? = 0, 

(5) 2t^ + rh^g = 0, 

(6) - Atl + 5tl + r^f°^ + rdrtl + ^4 = 0, 

(7) 2rH'i + 34 - rdrt'^ + 4tl + yr^t^ = 0, 

(8) ^rHl- 44 -rdrtl -34^0, 

(9) t%+^4 + Atl, + ^drt^T = 0, 

(10) 24 + ^rhl + 4tl-^dr4^0, 

(11) r{t'i + t^ -2rtl + drt% = 0, 

(12) ^rtl - rt? + 6rtl + 2dr4 + T^'^drtl = 0, 
5 ' 5 

(13) -drtl - -rdrtl - —tl - - Qtl = 0, 
r 5 5 

(14) ml + - rh^ + 2rdrtl = 0, 

(15) tl + + rHl + - 2^2 + rdrtl^ = 0. 

5 
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Those equations that only involve t°,tf and tf are given below: 

(1) ti + ltl-j4 + ti, = o, 

(2) rdr4 - 8^1 - + ^4 + 2^10 = 0, 

(3) dr4 + 2rtt^ + r f^tj - hlj ^ 0, 

(4) 4:r4 + drtl + rt"^ = 0, 

(5) 44 + rHl + ^4 = 0, 

(6) 64o - 4:4 + 5*9 + ^Vii + rdr4 + ^t? = 0, 

(7) 2rHl - + 34 - rdr4 + ^4 + ^r^4 = 0> 

2 5 

(8) ^rHt-d.4-Tdr4-34 + ]-4 = i), (E-2) 

(9) + ^^9 + 6^11 + -dr4 = 0, 

(10) 24 + ^r^4 ~ ^4 + 84 - ^dr4 = 0, 

(11) r{tl + 4]-Ar4 + dr4o = ^, 

(12) ^r4 - rtj + Sr4 + + '2dr4 + ^r^dr4 = 0, 

(13) -dr4 - -rdr4 - —4 -4- ^4 = o> 

7* 5 5 

(14) 104 + 84 - rhl + + 2rdr4 = 0, 

(15) 4 + \r^4 + rHl + - 44 + ra^t^i = 0. 

5 
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